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(57) Abstract 

A system (20) of the present invention includes light source instrumentation (30) to selectively illuminate a light scattering medium 
including a luminophore and detection instrumentation (50) to detect multiply scattered light output from the medium in response to 
illumination by the light source instrumentation (30). A processor (70) is operatively coupled to the detection instrumentation (50) to 
determine a first optical characterization of the medium from a first multiply scattered light output of a first illumination light wavelength 
and a second optical characterization of the medium from a second multiply scattered light output of a second illumination light wavelength 
different than the first illumination light wavelength. The processor (70) is operable to calculate lifetime of the luminophore from the 
first optical characterization, the second optical characterization, and a multiply scattered emission of the luminophore from the medium in 
response to excitation. 
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WO 00/2241 4 PCT/US99/23709 
CHARACTERIZATION OF LUMINESCENCE IN A SCATTERING MEDIUM 



BACKGROUND 

20 The present invention relates to spectroscopic techniques involving luminescence, 

and more particularly, but not exclusively relates to the determination of lifetime of a 
luminophore in a light scattering medium. 

There has been significant development of fluorescent and phosphorescent dyes or 
probes with decay kinetics dependent upon the presence or concentration of an analyte or 

25 metabolite. Accordingly, lifetime of these dyes can be measured to detect corresponding 
analyte(s) and/or metabolite(s) concentration. Fluorescent probes in the near-infrared 
range appear particularly promising for in vivo biomedical diagnostic techniques that 
involve external, noninvasive measurements or minimally invasive, endoscopic 
measurements of emitted light. 
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Unfortunately, quantitative lifetime measurements for such probes are often 
difficult to obtain in the light scattering environment typically encountered with in vivo 
diagnostics. Light scattering also hampers other applications of luminophore probes both 
inside and outside the biomedical field. Consequently, lifetime measurements are usually 
restricted to dilute, nonscattering solutions. Notably, even equipment used in this manner, 
such as a curvette to contain the dilute solution, tends to scatter light to some degree 
introducing an attendant inaccuracy. 

Moreover, current lifetime measurement approaches have other limits - especially 
for fluorophore probes. For example, deconvolution of instrument function often hampers 
accurate time-domain measurement of lifetimes. In another example, frequency-domain 
approaches generally require a reference fluorophore with known lifetime characteristics 
in the environment of interest and at the appropriate excitation and emission wavelengths. 
Thus, there is a need for further contributions that address these limits and/or other 
drawbacks confronting this technology. 
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SUMMARY OF INVENTION 

Accordingly, one form of the present invention is a unique technique to evaluate a 
medium including a luminophore. Other forms include unique systems and methods to 
measure optical properties of a luminophore in a light scattering medium. 

In another form, a light scattering medium includes a luminiphore that is exposed 
to a number of different wavelengths. Multiply scattered light from exposure to these 
wavelengths is measured, and one or more optical characteristics of the medium are 
determined relative to the different wavelengths. 

In still another form, an optical characteristic of a luminophore in a light scattering 
medium is determined. For this form, the medium is illuminated by light at a first 
wavelength corresponding to an emission wavelength of the luminophore and at a second 
wavelength corresponding to an excitation wavelength of the luminophore. Multiply 
scattered light in response to illumination by the different wavelengths is detected to 
optically characterize the medium and/or luminophore. 

A further form of the present invention is a technique to determine the lifetime of a 
fluorophore. This technique includes both a method and instrumentation directed to 
fluorescence lifetime measurements. Preferably, this technique does not utilize a reference 
fluorophore and is not adversely impacted by any light scattering and absorption that 
might occur in a scattering medium. 

Still a further form of the present invention includes a frequency-domain approach 
to measuring fluorescence lifetime of one or more fluorophores. This approach may 
include exciting a fluorophore in a light scattering sample with a modulated excitation 
light and detecting phase shift information. The fluorescence lifetime may be determined 
from the phase shift information through relationships characterizing light scattering by 
the sample. 

In yet a further form, a light scattering sample containing a fluorophore is exposed 
to an intensity modulated light at a first wavelength selected to cause the fluorophore to 
fluoresce at a second wavelength of light. Scattered light at the first wavelength is 
detected and scattered light at the second wavelength is detected. Optical properties are 
determined from the detected scattered light to provide fluorescence lifetime based on 
relationships that characterize photon migration in a light scattering medium. This form 
may include characterizing the detected scattered light in lerms of phase shift in the 
frequency domain. 
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In an additional form* a system of the present invention includes an intensity 
modulated excitation light source configured to deliver an excitation light of a first 
wavelength to a light scattering substance containing a fluorophore of interest. The 
fluorophore responds to the excitation light to provide a light emission at a second 
wavelength. Scattered light is detected at two locations spaced apart from one another. A 
processor gathers information corresponding to the detected scattered light and processes 
this information to determine fluorescence lifetime by applying relationships that 
characterize photon migration of scattered light. For this form, a first detector may be 
used for detection of light at a first one of the locations that includes a first optical fiber 
coupled to a first sensor. Also a second detector may be used for detection of light at a 
second one of the locations that includes a second optical fiber coupled to a second sensor. 
The first and second detectors may also include corresponding first and second optical 
filters to selectively detect the first wavelength of light with the first sensor and the second 
wavelength of light with the second sensor. The coupling arrangement of the first and 
second fibers to the first and second sensors may be interchanged to obtain comparative 
measurements for the minimization of equipment inaccuracies. 

Further forms, embodiments, objects, features, aspects, advantages, and benefits of 
the present invention shall become apparent from the detailed description and drawings of 
the present application. 



WO 00/22414 



PCT/US99/23709 



BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is a schematic view of a system according to one embodiment of the present 
invention. 

Figs. 2 A and 2B depict a flow chart illustrating one process that may be performed 
5 with the system shown in Fig. 1 . 

Fig. 3 is a schematic view of a system of another embodiment of the present 
invention. 

Fig. 4 depicts normalized excitation and emission spectra for DTTCI (top) and 
ICG (bottom) in a comparative format pertinent to experimental examples. 
10 Fig. 5 plots phase shift versus modulation frequency for two different 

concentrations of DTTCI in a comparative format pertinent to experimental examples. 

Fig. 6 plots phase shift versus modulation frequency for two different 
concentrations of ICG in a comparative format pertinent to experimental examples. 

Fig. 7 plots lifetime versus modulation frequency for DTTCI and ICG 
15 corresponding to the plots of Figs. 5 and 6. 

Fig. 8 plots relative modulation attenuation versus modulation frequency for two 
different concentrations of DTTCI in a comparative format pertinent to experimental 
examples. 

Fig. 9 plots relative modulation attenuation versus modulation frequency for two 
20 different concentrations of ICG in a comparative format pertinent to experimental 
examples. 

Fig. 10 plots lifetime versus modulation frequency for DTTCI and ICG 
corresponding to the plots of Figs. 8 and 9. 



5 
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DESCRIPTION OF PREFERRED EMBODIMENTS 
For the purpose of promoting an understanding of the principles of the invention, 
reference will now be made to the embodiments illustrated in the drawings and specific 
language will be used to describe the same. It will nevertheless be understood that no 
5 limitation of the scope of the invention is thereby intended. Any alterations and further 
modifications in the described embodiments, and any further applications of the principles 
of the invention as described herein are contemplated as would normally occur to one 
skilled in the art to which the invention relates. 

As used herein, "lifetime" refers to the mean survival time of an activated 
10 luminophore or the mean time between the absorption of an excitation photon and 

emission of a photon. Further, as used herein "multiply scattered light" refers to light that 
travels at least five (5) times the mean isotropic scattering length [(l-g)^]" 1 ; where g is 

the mean cosine of angular scatter and |i s is the scattering coefficient of the medium. 

Fig. 1 illustrates evaluation system 20 of one embodiment of the present invention. 

15 System 20 includes light source instrumentation 30, container 40, detection 

instrumentation 50, and processor 70. Light source instrumentation 30 includes 
modulation signal generator 22 having a range of selectable output Radio Frequencies 
(RF). Generator 22 drives two monochromatic light sources in the form of laser diodes 
24, 26. Laser diodes 24, 26 are intensity modulated at a selected RF frequency (co) input 

20 from generator 22 to provide light at an excitation wavelength (XJ and emission 

wavelength (Xm), respectively. Instrumentation 30 also includes kinematic mirror 32 that 
is operable to successively select between the two laser beams at wavelengths X x , Xm 
output by the respective laser diodes 24, 26. Light from mirror 32 encounters 
continuously variable neutral density filter wheel 34 of instrumentation 30 to selectively 

25 adjust intensity. Lens assembly 36 of instrumentation 30 collects the light output by 
neutral density filter wheel 34 for input to optical source fiber 38. 

Source fiber 38 enters container 40. Container 40 holds a light scattering medium 
M including a selected amount of a luminophore as a constituent. Source fiber 38 
discharges light at site 39 of medium M corresponding to the origin position r<>(r=0). 

30 Detection instrumentation 50 includes optical detector fibers 44a, 44b having respective 
light input sites 42a, 42b in medium M. Sites 42a, 42b are spaced apart from each other 

6 
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and correspond to radial distances r t and r 2 relative to r 0 ; where this spacing is represented 
by At (Ar = r? - T\). 

Interchangeable connectors 46a, 46b couple fibers 44a,44b to detector 48 of 
detection instrumentation 50. Detector 48 includes two optic channels 50a, 50b coupled 

5 by connectors 46a, 46b to receive light from detector fibers 44a, 44b, respectively. Each 
channel 50a, 50b has a corresponding interchangeable/removable interference filter (IF) 
54a, 54b; adjustable neutral density filter (NDF) 56a, 56b; and light sensor 58a, 58b. 

Light sensors 58a, 58b are coupled to processor 70 to provide one or more output 
signals corresponding to light detected from medium M. Sensors 58a, 58b are arranged in 

10 a standard heterodyne configuration and may be of any form such as Photomultiplier 
Tubes (PMTs), photodiodes, or image intensified Charge Coupled Devices (CCDs) to 
name just a few. For the heterodyne configuration, RF signal generator 60 is phase locked 
to generator 22 at a slightly different frequency co + Aw as represented by coupling 62; 
where Aco is the frequency difference. The output of generator 60 is amplified by RF 

15 amplifier 64 and mixed at sensors 58a, 58b to provide a corresponding differential output 
signal from which information corresponding to phase and modulation magnitude of the 
detected light can be determined with processor 70. Accordingly, processor 70 is also 
operatively coupled to generators 22, 60. 

Processor 70 includes a port for insertion and removal of a portable memory device 

20 such as an electromagnetically or optically encoded disk, cartridge, or tape. In addition to 
being coupled to light source instrumentation 30 and detection instrumentation 50, 
processor 70 is also operatively coupled to visual display 76 and one or more Input/Output 
(I/O) devices 78, including for example a keyboard, mouse, light pen, acoustic 
loudspeakers, microphone, and/or printer just to name a few. Processor 70 may be 

25 comprised of one or more components configured as a single unit, or when of a multi- 
component form, processor 70 may have one or more components remotely located 
relative to the others, or otherwise have its components distributed throughout system 20. 
Processor 70 may be programmable, a state logic machine or other type of dedicated 
hardware, or a hybrid combination of programmable and dedicated hardware. One or 

30 more components of processor 70 may be of the electronic variety including digital 

circuitry, analog circuitry, or both. As an addition or alternative to electronic circuitry, 
processor 70 may include one or more optical elements. 
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Processor 70 includes an integrated and/or remote storage capability in the form of 
one or more types of memory. By way of nonlimiting example, this memory may include 
one or more of the solid-state, magnetic, and/or optical memory types. Such memory 
types may include Random Access Memory (RAM), Sequential Accessible Memory 
5 (SAM) (such as the First-In, First-Out (FIFO) variety, or the Last-In, First-In LIFO 

variety), Programmable Read Only Memory (PROM), Electrically Programmable Read 
Only Memory (EPROM), flash memory or Electrically Erasable Programmable Read Only 
Memory (EEPROM); an optical disc memory (such as a CD ROM); a magnetically 
encoded hard disc, floppy disc, tape, or cartridge; another variety of storage device as 

10 would occur to those skilled in the art, or a combination of any of these types. 

Furthermore, the memory may be volatile, nonvolatile, or a hybrid combination of volatile 
and nonvolatile varieties. Also, memory may be permanently installed, in a portable form 
that may be readily removed and reinstalled, or a combination of these types. 

In one embodiment including electronic circuitry, processor 70 is of a standard 

15 personal computer configuration with a common solid-state digital integrated processing 
unit operatively coupled to solid-state memory. For this personal computer embodiment, 
appropriate interfaces are installed to facilitate control of generators 22, 60 and receipt of 
data from detector 48. The memory of this embodiment contains programming to be 
executed by the processing unit, and is arranged for reading and writing of data in 

20 accordance with one or more routines executed by processor 70. Besides memory, 
processor 70 may include any oscillators, control clocks, interfaces, signal 
compensators/conditioners, filters, limiters, Analog-to-Digital (A/D) converters, Digital- 
to-Analog (D/A) converters, communication ports, or other types of circuits as would 
occur to those skilled in the art to implement the present invention. 

25 Processor 70 is configured to execute one or more routines to perform selected 

calculations with data received from detector 48 for lifetime evaluation process 120. 
Referring additionally to the flow chart of Figs. 2A and 2B, evaluation process 120 is 
illustrated. Evaluation process 120 utilizes Frequency Domain Photon Migration (FDPM) 
techniques to extract the characteristic optical properties of medium M. These properties 

30 are obtained for photons of two wavelengths: the wavelength used to optically excite a 
selected luminophore and the wavelength of the luminescence resulting from this 
excitation. Lifetime measurements are determined from these wavelength-dependent 
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10 



15 



20 



25 



characterizations of the medium and comparative luminescence measurements, as 
described hereinafter. 

In stage 122 of process 120, a luminiphore probe with known excitation 
wavelength X x and emission wavelength X m is selected for evaluation and placed in light 
scattering medium M. This light scattering medium M may include, for example, living 
biologic tissue for which metabolites/analytes are being interrogated in terms of lifetime of 
the selected luminophore probe. In other examples, the light scattering medium M may be 
a cell culture, flow cytometry stream, a chemical reaction medium or other light scattering 
environment as would occur to those skilled in the art. In one alternative, endogenous 
luminophores in medium M are interrogated without introduction of an exogenous probe 
in stage 122. 

In stage 130, intensity modulated light at frequency-co with the excitation 
wavelength X x of the designated luminophore is provided by light source instrumentation 
30 to source site 39 of medium M. Accordingly, light from laser diode 24 is provided 
from light source instrumentation 30 through fiber 38 to site 39. The illuminating light is 
subsequently scattered and/or absorbed by medium M. As multiply scattered light sourced 
from site 39 reaches sites 42a, 42b; it can be sensed with detection instrumentation 50 and 
quantitized in terms of frequency domain parameters of relative phase shift and/or 
modulation magnitude. Filters of detector 48 are removed and/or adjusted to facilitate 
detection of the excitation light wavelength with sensors 58a, 58b for this stage. 

Photon transport in a light scattering medium M may be modeled as a diffusive 
process. In the frequency domain, the photon density U(r,co) in a homogeneous medium 
at a vector position r can be related to optical properties of the medium by the diffusion 
equation ( 1 ) as follows: 



where: D = [3(n' s + ^a)]* 1 is the diffusion coefficient of the medium, c is the speed of light 
in the medium, q(r,co) describes properties of the light source, and co is the modulation 
angular frequency of the light source (generator 22). The output signal from detection 
instrumentation 50, V(r,co), depends on the complex responsivity R(X,co) expressed as: 
V(r,co) = R(X,co)U(r,co). The magnitude of R(X,co) represents gain and conversion 



-cDV 2 U(r,co) + (en, + /co)U(r,co) = q(r,o); 



(1) 
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efficiency of the photon transport and the angle represents the phase delay of the 
modulated light. 

For the diffusion equation model, the absorption coefficient ji a and the isop tropic 
scattering coefficient |T S characterize pertinent light scattering properties of the medium at 

5 a given light wavelength. The isotropic scattering coefficient is related to the scattering 
coefficient m by: ja\ = (l-g)n$; where g is the mean cosine of angular scatter. As used 
herein, the subscripts tfc x" and "m" are used to designate various optical parameters 
specific to the excitation and emission wavelengths, respectively. It has been found that 
lifetime measurements in a multiply scattering medium may be performed by accounting 

10 for the optical characteristics of the medium at two light wavelengths commonly 

associated with the luminophore: (a) the excitation wavelength X x and (b) the emission 
wavelength X m . This characterization can be in terms of the absorption and isotropic 
scattering coefficients for each of the two wavelengths: (a) ja ax , i±\ % for X x and (b) > 
Ji'sm for X m . 

15 Accordingly, in stage 130 medium M is characterized at X x beginning with 

operation 132. In operation 132, medium M is exposed to modulated light of wavelength 
X x at site 39 as sourced from laser diode 24 of light source instrumentation 30. This source 
light can be modeled as an isotropic point source of the form: q x (r,co) = P x (a>)5(r); where 
P x (co) is a complex number that represents the source magnitude and phase and 5(r) is the 

20 Dirac delta function. Furthermore, by modeling the light behavior in terms of the 

diffusion equation with infinite boundary conditions, the solution to the diffusion equation 
takes the form of spherical photon density waves described by the following expression 
(2): 



25 



PA") 
4 7T cD x r 



(2) 



where: U x (r,co) is the frequency domain excitation photon density in the medium M and 
30 the complex wave vector /: x (co) is given by expression (3) as follows: 
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(3) 



On substitution of expression (3) into expression (2), the excitation photon density can be 
modeled as follows in expression (4): 



where: 

\l/2 



(4) 



(3 x {u>) = a x (oj) cos^ tan" 1 ^^— 
7x(w) = a c (u;) sin ^ tan " l ^~"j) 



The photon density is related to the observed modulation phase G(r,co) by 
expression (5a) as follows: 



♦ at \ ImL7(r ' u) (5a) 



The modulation M(r,co) of the photon density waves a distance r away from and 
normalized to unity at the source (r = 0) is related to photon density by expression (5b) as 
follows: 



M(r,w) = 



Tie* U(r,u>) + lm* U(r,u) (5b) 
\] 1le*U{r,Q)+Im*U(r,Q) ' 



For injection of source light at the excitation wavelength into the medium M in operation 
132, the substitution of expression (4) into equations (5a) and (5b) results in expressions 

it 
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(6a) and (6b), respectively, as follows: 

G x (r,co) = -/ x (co)r (6a) 
M t (r,u,) = ^e^<°)-^. (6b) 

The relative modulation phase and magnitude are observed with detection 
10 instrumentation 30 at the two radial distances n and r 2 ; where Ar = r 2 - n t T2 > n. The 
resultant phase difference is expressed as: A9 x (Ar,co) = y x (co)Ar; and the resultant 
modulation magnitude may be as follows in expression (7): 



M x (r U Lj) 



where the actual modulation information observed with detection instrumentation 50 is as 
follows: 

20 

m r (r,u/) = M = (t,uj) m sx m <ix (r,cj) (8) 



25 and is related to M x (r,co) by the modulation of source m sx and the modulation response 

nvj x (r,co) of detection instrumentation 50. The ratio of the observed modulation signals is 
provided by expression (9) as follows: 



This ratio is generally independent of source modulation. 

12 
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In operation 132, data is collected and stored in processor 70 corresponding to 
detected multiply scattered light output from medium M. It should be appreciated that 
only one of relative phase and modulation attenuation information needs to be observed to 
provide the desired characterization of the medium M at the excitation wavelength. In 
5 order to calculate both the absorption and isotropic scattering coefficients ja^ , ^ sx with 
expression (4) from relative phase or magnitude observations, a regression or other 
iterative estimation can be employed. To facilitate this calculation, operation 132 detects 
multiply scattered light outputs for a number of different RF modulation frequencies. 
While relative measures of phase or modulation attenuation are indicative of 

10 photon-migration through medium M, differences in the response function of channels 50a 
and 50b of detector 48 may result in inaccuracy that reaches an undesirable level for some 
applications. It has been found that instrumentation response function effects may be 
reduced by switching from the illustrated (first) configuration with sites 42a, 42b 
respectively coupled to channels 50a, 50b to a second configuration with sites 42a, 42b 

15 coupled to channels 50b, 50a; and repeating the phase and/or modulation magnitude 
measurements of operation 132 in operation 134 with the second configuration. The 
measurements of operation 134 are conducted with the same relative spacing between sites 
39, 42a, and 42b. Connectors 46a, 46b provide a convenient way to manually perform this 
reconfiguration. In an alternative embodiment, an optical multiplexer may be utilized to 

20 automatically accomplish this reconfiguration under the control of processor 70. 

In an example based on relative phase measurements, the effective response 
function of the detection instrumentation is designated 8j nstr . Letting AOxforO = 9 x (r 2 ) - 
9x(i"i) + Oinstr represent measurements during operation 132 and A9 x (rir2) = 9 x (ri) - Gxfo) + 
Qinstr represent measurements during operation 134, the instrument effects can be removed 

25 by taking one half of the difference between these two relative measurements. The desired 
phase difference, A9 X , is then expressed as: A9 X = 0.5[A9 x (r 2 ri) - A9 x (rir 2 )]. 

Measurements gathered with detection instrumentation 50 in operations 132, 134 
of stage 130 for a desired number of different modulation frequencies are recorded with 
processor 70. Stage 140 is next encountered to obtain an optical characterization of 

30 medium M at the emission wavelength X m in terms of absorption and isotropic scattering 
coefficients ^i am , jT sm . In operation 142, emission wavelength light is provided to source 
site 39 with laser diode 26 of light source instrumentation 30, instead of excitation 
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wavelength light from laser diode 24. Filtering of detector 48 is removed/adjusted to 
facilitate detection of multiply scattered light output from medium M at the emission 
wavelength X m and measurement of relative modulation phase or magnitude over a 
selected range of modulation frequencies in the manner described for stage 130. The 
emission photon density is characterized in stage 140 in accordance with expressions (2)- 
(9) substituting the emission wavelength for the excitation wavelength (and 
correspondingly substituting subscript "x" with "nT in these expressions). In operation 
144, the measurements are repeated with coupling of site 42a being switched to sensor 58b 
of channel 50b and coupling of site 42b being switched to sensor 58a of channel 50a. The 
measurements of operations 142, 144 are recorded with processor 70 to determine 
absorption and isotropic scattering coefficients |i a m» n'sm- 

In stage 150, luminescence data is gathered. Excitation wavelength light is 
provided to site 39 of medium M from laser diode 24 of light source instrumentation 30 in 
operation 152 of stage 150 and frequency domain measurements are gathered. For stage 
150, detection instrumentation 50 is configured with filtering of channel 50a 
adjusted/removed to detect the excitation wavelength with sensor 58a and filtering of 
channel 50b adjusted/removed to detect the emission wavelength with sensor 58b. For a 
fluorescent type of luminophore, fluorescence photon density may be modeled in 
accordance with expression (10) as follows ( where the subscript "f denotes fluorescence 
optical parameters): 

rr , ^ P M (e-^'-e-^ \ f l+iun\ (l0) 

The quantum efficiency of the fluorophore is denoted by #and n af describes absorption of 
excitation light due to fluorescence. Fluorescence decay is assumed to be of the 
monoexponential type with lifetime x; however, the principles of the present invention can 
be applied to multiexponential decays using techniques known to those skilled in the art. 
Using expression (3), the fluorescence photon density of expression (10) can be written as 
presented in the following expression (11): 
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6^((d) (11) 
U f (r, id) = f — -p- {(>(r, ~ K(r, co)u>t] 

where: 

&(>. w)i -r c(r, u>)p(u>) 

K(r,a,) = F+IpM] 5 • 

6(r f a)) - expl-p.Cuiirlcos^JuiJr] 

- cxpt-p m (w)r]cos[7 m (ujir] f 
£(r, to) = exp[-p x (u*)rlsin[-y,(o>)r] 
10 - exp[-p m (aj)r]sin[7mMr], 

6 D m D/ 
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il> / 1 1 \ 



15 The fluorescence photon density is related to the observed fluorescence modulation phase 
6Kr,co) by expression (5a). Substituting expression (11) into expression (5a), expression 
(12) results as follows: 

K(r f o>) + Mr, u>Uut 

tan 6,(r, w) - — ■ — , (1^) 

U/(r, <d) — K(r t oj)cjt 



and the fluorescence decay lifetime can be written per expression (13) as follows: 



1 tan Bf\r. to) - T)(r, id) 
<d ^fr, u))tan 8,(r, id) + 1 * 

25 

where r\{r 9 to) = K(r, ca)A|/(r,a>) and determinations of i|/(r,co) and K(r, co) are based on 
measured phase shifts A9 x (Ar,co) and A0(<Ar,co). 

Alternatively, the fluorescence photon density is related to fluorescence 
modulation MKi\cd) by expression (5b). Substituting expression (11) into expression (5a) 
30 and noting that *(r, 0) = 0 yields expressions (14) and (1 5) as follows: 

c{r 2 ^) (14) 



p(o) ^ 
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where: 



e(r 2 ,u;) = 



^(r 2 ,0) 



(15) 



5 



10 



is based upon the optical coefficients of the sample. Comparable expressions can be 
derived for other types of luminescence and different boundary conditions. 

After measurements of relative modulation phase or magnitude are made in 



operation 152, the measurements are repeated in operation 154 with sites 42a, 42b being 
switched relative to channels 50a, 50b as described in connection with operation 134 and 
144, As part of the reconfiguration for stage 154, the filtering of channels 50a, 50b is also 
adjusted to facilitate detection of the excitation wavelength with sensor 58b and the 
15 emission wavelength with sensor 58a. Measurements of operations 152, 154 are stored 
with processor 70. 

In calculation stage 160 (Fig. 2B) of evaluation process 120, processor 70 performs 
calculations with the data collected during stages 130, 140, 150. In operation 162, 
processor 70 performs regression analysis of measurements from excitation 
20 characterization stage 1 30 to determine the absorption and isotropic scattering coefficients 
jiax , |i' sx at the excitation wavelength. In operation 164, processor 70 performs regression 
analysis of measurements from emission characterization stage 140 to determine the 
absorption and isotropic scattering coefficients n am , lAmat the emission wavelength. 



25 obtained from luminescence characterization stage 1 50 and the coefficients obtained from 
operations 162, 164. For calculations of lifetime based on relative phase measurements in 
stage 150, the measured phase shift AG t <Ar,co) = Qfa ,co) - 9 x (ri ,co) is adjusted by the 
addition of phase change 0 xc (ri associated with the excitation propagation from site 39 
to site 42a at n. This phase value G xc (ri ,co) can be calculated from expression (6a) and the 

30 optical characteristics determined in operation 162 from measurements of excitation 
wavelength characterization stage 130. The resultant modulation phase 9f(r2 ,co) - 



In operation 166, processor 70 calculates lifetime based on the measurements 
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A9*(Ar,co) + 0 xc (ri *<o) is the phase at the source (r = 0) and can be used in expression (13) 
to determine lifetime t. 

For calculation of lifetime based on modulation attenuation, the modulation 
magnitude at radial distance r> is referenced relative to the excitation modulation 
magnitude at radial distance n. The ratio of expression (14) to expression (6b) yields 
expression (16) and (17) as follows: 



M f (r 2 ,u>) 



M y (Ar,u;) - - c r l (..(o)- fcl w)] ^ (1 + (u , r )2) (16) 



c(r 2 ,u>) 



from lifetime x is derived: 



15 



t(r 2 ,u) 



- 1 



(17) 



20 



The modulation magnitude detected with detection instrumentation 50 is given by 
expression (18) as follows: 



m/{ r 2,w) = M/{r 2i u) m 5T m t i m (r 2 ,u;) 



(18) 



which is related to modulation Mf(r,co) by modulation m sx of the source at the excitation 
wavelength, and the modulation mdmfo of the detection instrumentation 50 at the 
25 emission wavelength. The ratio of expression (18) and the excitation modulation of 
expression (8) is provided by expression (19) as follows: 



30 



x rnj(r 2 ,Lj) M;(r 2 ,uQ m dm (r 7i u>) 

m/(Ar, u;) = = -77-7 : 

;v m x (r u u> M x (r u uj m dx {r u uj 



(19) 



From the ratio of expression (19), the following expression (20): 
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MA&r^) = rn f (Ar^) H^h^l (20) 

rn dm (r 2 ,u> 



is found to depend upon the ratio of the detection instrumentation 50 response functions at 
excitation and emission wavelengths. This information may be expressed as a single ratio 
of modulations of the sources; where given the optical coefficients of the medium M from 
operations 162 and 164, excitation M xc (r,a>) and emission Mmc(r,co) can be determined 
from expression (7). The product of source modulation and detection instrument response 
function can be determined from expression (8) at the excitation: 



m r (ri) 

m« m^u) = Mcc{ri u) (2D 



and emission: 

/ ^ ( 22 ) 

m, m m dm (r 2 ,u,)- 

wavelengths. The resulting ratio of detection instrument response functions is given by 
expression (23) as follows: 

ra<£x(ri,u;) ^ m sm m g (r 1 ,u;) M mc (r 2 ,u;) ^ 

m dm (r 2 ,u) m ST m m (r 2 ,u;) M xc (r x ,u) 



which depends upon a comparison between measured and calculated modulation 
information and a constant ratio of modulations of the sources at emission and excitation 
wavelengths. Varying the source modulation ratio, the lifetime measured at multiple 
modulation frequencies can be regressed to obtain a unique source modulation ratio 
associated with a vanishing slope, i.e. minimized x* of the lifetime distribution, and to 
obtain the resulting lifetime x. After lifetime is determined in stage 160 from relative 
phase or modulation measurements, it is output in stage 170, concluding process 120. 
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It should be understood that only one possible sequence of the measurement stages 
130, 140, 150 is illustrated. Indeed, the various measurements may be performed in many 
other sequences with respect to selected light wavelengths, modulation frequencies, and 
coupling configuration of sites 42a. 42b relative to channels 50a, 50b. Also, the 
5 calculations of stage 160 may be performed any time by processor 70 in relation to the 
measurement stages 130, 140, 150 to the extent measurement data has been provided to 
processor 70. 

Ln other embodiments, the arrangement of system 20 may differ. For example, 
light source instrumentation may include a source of excitation and/or emission 

10 wavelength light such as a different type of laser, lamp, or other device as would occur to 
those skilled in the art. In another example, detection instrumentation 50 may include an 
optical multiplexer controlled by processor 70 to provide for the switching of optical 
couplings relative to detector 48. In one alternative embodiment, the measurements may 
be based on two spaced apart light source sites 39 in medium M with just one detection 

15 site 42a or 42b and corresponding sensor 58a or 58b. In one form of this arrangement, 
light source instrumentation may be configured with another set of the components 
described for instrumentation 30 to provide two light source channels; while detection 
instrumentation 50 could be modified to include only a single channel 50a or 50b. The 
lifetime calculations may be readily adapted to this dual source using techniques known to 

20 those skilled in the art. Also, instrument function that might lead to calculation 

discrepancies relative to the tw o source sites can be addressed by switching the equipment 
supplying light to each source; thereby providing the first and second measurement 
configurations described in connection with stages 130, 140, 150. 

Fig. 3 depicts medical diagnostic system 220 of another embodiment of the present 

25 invention. System 220 includes diagnostic instrument 225 to evaluate a patient's medical 
condition based on lifetime readings of an endogenous, exogenous, or immobilized 
exogenous fluorophore in a patient's tissue 280. Instrument 225 includes light source 
apparatus 230, detection apparatus 250, and processor 270 that are operationally 
configured like light source instrumentation 30, detection instrumentation 50, and 

30 processor 70, respectively, but are packaged in a manner convenient for in vivo 
interrogation of tissue 280. 

The fluorophore in tissue 280 is typically an exogenous probe introduced into 
tissue 280 that has a known excitation and emission wavelength to which the optics of 
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instrument 225 are matched. Such exogenous probes may be immobilized in a 
subcutaneous implant, by oral ingestion, or any other means that would occur to those 
skilled in the art. Alternatively or additionally, the present application may be applied to 
monitor endogenous fluorophores for which excitation/emission wavelengths are known. 

Light source apparatus 230 includes one or more sources to provide light to tissue 
280 at excitation and emission wavelengths for the selected fluorophore in accordance 
with evaluation process 120. The source light is provided to site 239 of tissue 280 as 
represented by arrow IL. The source light is multiply scattered by tissue 280 and received 
by detection apparatus at sites 242a, 242b as represented by arrows MSI, MS2. 
Instrument 225 is configured to maintain a generally constant relative spacing between 
sites 239, 242a, and 242b, corresponding to r 0 , n, and r 2 , respectively. Processor 270 may 
be configured to control selection of the appropriate source light wavelength and 
modulation frequencies for light source apparatus 230 and the corresponding filter 
configuration for detection apparatus 250 to automatically perform the stages and 
operations of process 120 with instrument 225. 

For external interrogations, instrument 225 may include a moveable hand-held 
wand or optical probe coupled to a base unit; where the wand is arranged for placement 
proximate to surface 282 of tissue 280. In this case, surface 282 may be the patient's skin, 
with instrument 225 being arranged to make percutaneous measurements. In other 
instances, instrument 225 may include an endoscope coupled to a base unit that is arranged 
to perform interrogations through a body lumen, cavity, or small surgical incision. For 
these instances, surface 282 can be the boundary of the body lumen or an organ selected 
for interrogation, to name just a few examples. Depending on the configuration of 
instrument 225, the calculations performed in stage 160 may be adjusted to account for 
finite boundary conditions to improve performance. In another embodiment, instrument 
225 is arranged to place sites 239, 242a, 242b in the tissue, more closely approximating 
infinite boundary conditions. Furthermore, in alternative embodiments, instrument 225 
may be completely contained in a portable device that is battery powered. 

Yet another embodiment includes interrogating a light scattering medium including 
an amount of a selected luminophore with light at an emission wavelength of the 
luminophore and sensing multiply scattered light at the emission light wavelength in 
response to this interrogation to provide a first optical characterization of the medium. 
The medium is also exposed to light at an excitation wavelength of the luminophore and a 

20 
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multiply scattered luminescence at the emission wavelength is detected in response to this 
exposure. A value corresponding to lifetime of the luminophore is determined from the 
first optical characterization and the luminescence. 

A further embodiment includes exposing a light scattering medium including an 
5 amount of a luminophore to a number of different light wavelengths selected relative to 
the luminophore. A plurality of optical characteristics of the medium are established by 
sensing multiply scattered light at each of the different wavelengths. A value 
corresponding to lifetime of the luminophore is determined from the optical characteristics 
and a multiply scattered emission at a first one of the light wavelengths caused by 
10 exposure of the medium to a second one of the light wavelengths. 

Still a further embodiment provides for an evaluation of a light scattering medium 
including an amount of a luminophore with a light interrogation system. This system 
includes a first device optically coupled to a first site of the medium and a second device 
optically coupled to a second site of the medium that is spaced apart from the first site. 
15 Also included are subjecting the medium to a first light wavelength selected relative to the 
luminophore, optically coupling the first device to the second site and the second device to 
the first site, and exposing the medium to the first light wavelength after coupling. A 
value corresponding to lifetime of the luminophore is determined from a first light output 
and a second light output. 
20 Another embodiment includes means for illuminating a light scattering medium 

including a luminophore; means for characterizing light scattering behavior of the medium 
for an excitation wavelength of the luminophore and an emission wavelength of the 
luminophore; and means for determining lifetime of the luminophore from the 
characterizing means and a multiply scattered light emission from the medium at the 
25 emission wavelength in response to illumination by light at the excitation wavelength. 

Also, because luminophores may have different activated states, they may have 
multiple lifetimes corresponding to these states. While the examples described herein are 
directed to the lifetime of luminophores that are excited to a single activated state to 
preserve clarity, in other embodiments the principles of the present invention can be 
30 applied to luminophores with multiple activated states and lifetimes using techniques 

known to those skilled in the art. Such multiple lifetime probes are also within the scope 
of the present invention. 



21 



WO 00/22414 



PCT/US99/23709 



EXPERIMENTAL EXAMPLES 
The present invention will be further described with reference to the following 
specific examples. These experiments and results are intended to be illustrative of the 
present invention and should not be considered limiting or restrictive with regard to the 
5 scope of the present invention. Further, any theory, mechanism of operation, proof, or 
finding stated herein is meant to further enhance understanding of the present invention 
and is not intended to make the present invention in any way dependent upon such theory, 
mechanism of operation, proof, or finding. 

Several experiments were conducted with a test setup comparable to system 20. 

10 FDPM measurements according to process 120 were conducted in a light scattering tissue- 
simulating phantom form of medium M. For these experiments the lifetimes of different 
micromolar concentrations of two luminophores of the fluorescent type where evaluated: 
(1) 3,3'-Diethylthiatricarbocyanine Iodide (DTTCI) and (2) Indocyanine Green (ICG or 
IR-125) (both from AR-COS Organics, Fair Lawn, N.J.). Container 40 for the phantom 

15 was a cylindrical acrylic vessel 1 1 .3 cm in diameter and 15 cm in height. The phantom 
was prepared as an aqueous (DUIF water; Fisher Scientific, Fair Lawn, N J.) intralipid 
(20%; Pharmacia & Upjohn Company, Clayton, N.C.) solution to mimic tissue-like 
scattering properties. A DTTCI stock solution was prepared with ethanol; and ICG was 
dissolved directly in water. Appropriate quantities of stock solutions were suspended in 

20 1 .0% Intralipid to yield the final fluorophore concentrations in the phantoms as listed in 
Table 1 that follows: 



Dye 


Cone. 










Mam 








fc,M] 


(nmj 


[1/cm] 


[1/cm] 


[nm] 


[1/cmj 


[1/cm] 


[ns] 


DTTCI 


0.5 


749 


0.054(S) 


9.2(12) 


828 


0.032(3) 


10.2(8) 


1.34(3) 




1.0 


749 


0.07(2) 


6.8(15) 


828 


0.031(3) 


S.4(7) 


1.34(4) 


ICG 


0.0625 


778 


0.039(8) 


9.5(13) 


S2S 


0.033(4) 


S.0(S) 


0.54(3) 




0.125 


77S 


0.05(1) 


7.9(20) 


828 


0.051(9) 


9.2(15) 


0.56(4) 
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Additional samples were prepared in water for spectral analysis. Figure 4 shows the 
concentrations of the analyzed dyes and presents the relevant excitation and emission 
spectra measured with a spectrofiuorometer (Fluorolog-2; SPEX Industries, Inc., Edison, 
N.J.). 

The lifetime of DTTCI dissolved in ethanol was previously reported as 1.33 ± 0.02 
ns at 790-nm excitation and 820-nm emission wavelengths. The lifetime of ICG was 
reported to be 0.57 ± 0.02 ns at 780-nm excitation and 830-nm emission wavelengths 
measured without scattering by use of time-domain techniques as well as by conventional 
frequency-domain fluorimetry referenced against the known lifetime of a well- 
characterized fluorophore. 

For the experimental examples, excitation wavelengths were provided by laser 
diodes of the 56 DFS series; Melles Griot, Boulder, Colorado with power/wavelengths of 
3 mW at 749nm for ICG excitation and 25 mW at 778nm for DTTCI excitation. A 30 
mW 830 nm laser diode from Melles Griot was used for emission wavelengths. Kinematic 
mirror 32 was a model 9891 from New Focus, Santa Clara and neutral density filter wheel 
34 was from Newport of Irvine. California. Source fiber 38 was a 1000-|im-core source 
fiber (3M FT SILICA/0.39-NA TECS multimode fiber type). Also, fibers 44a, 44b were 
provided by two additional 1000-^m-core fibers of identical length placed in medium M at 
equal depth and parallel to source fiber 38. Detector fibers 44a, 44b were positioned 
within the sample at distances r : = 1 .0 cm and r 2 = 2.5 cm away from source fiber 38. The 
radial separation distance (Ar = 1.5cm) of the two detection fibers remained unchanged 
during all measurements. The light sensed by the fibers was delivered to sensors 58a, 58b 
in the form of gain-modulated photomultiplier tube (PMT) detectors (Model H6573; 
Hammamatsu, Bridgewater, N.J.). The detector fibers 44a, 44b were terminated with 
connectors 46a, 46b in the form of SMA fiber-optic connectors that could be interchanged 
conveniently to perform the reconfigurations of stages 130, 140, 150. 

Detector 48 included neutral-density filters 56a, 56b from CVI Laser Corporation, 
Albuquerque, N.M. and interference filters 54a, 54b of a narrow-bandpass form (lOnm 
FWHM; CVI Laser Corporation). Also included were lens assemblies that focused the 
collected light onto the PMTs. The gain settings of the PMTs remained unchanged during 
all measurements. The three neutral-density filters in the setup aided in maintaining 
constant dc levels of the detected PMT signals. The filters 34, 54a, 54b were used to adapt 
the setup to different output power levels of the two source laser diodes 24, 26 and to 
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facilitate acquisition of light intensity signals over a large dynamic range. The PMTs were 
gain modulated at the modulation frequency of the laser diodes plus an offset frequency of 
100 Hz (Aco). The heterodyned PMT signals were then digitized, Fourier transformed, and 
analyzed for phase shift and amplitude attenuation by use of Lab VIEW software (National 
5 Instruments Corporation, Austin, Tex.) executed by a personal computer form of processor 



The experimental results illustrated in Figs. 5-7 are based on relative phase 
measurements in stages 130, 140, 150. To obtain these results the stages were performed 
in a different order than illustrated in the flow chart of Figs. 2 A and 2B. Namely, Stage 

10 150 was performed first, followed by stage 130, and then stage 140. The measurement 
protocol for the lifetime determination comprises the three stages 130, 140, 150 for each 
concentration; where each stage includes measurements for two different configurations 
over a range of modulation frequencies from about 60 to about 140 MegaHertz (MHz). 
The kinematic mirror placed in the experimental setup expedited execution of the three 

15 stages of each measurement. Optical parameters of the sample at excitation and emission 
wavelengths were obtained from two-parameter least-squares fits of expression (5) to the 
excitation and the emission data, respectively. These optical coefficients together with the 
fluorescence phase shift then permitted the deduction of the lifetime of the fluorescent dye 
from expression (13). 

20 The top graph (a) of Fig. 5 provides experimentally determined plots of phase shift 

versus modulation frequency for a 0.5fiM concentration of DTTCI with the excitation 
characterization measurements of stage 130, the emission characterization measurements 
of stage 140, and the fluorescence characterization measurements of stage 150 shown in 
ascending order. The bottom graph (b) of Fig. 5 provides plots of the same type and in 

25 the same order as top graph (a) for a 1.0|xM concentration of DTTCI. The top graph (a) 
and the bottom graph (b) of Fig. 6 provide excitation, emission, and fluorescence 
characterization plots in the same order as the graphs of Fig. 5 for 0.0625|iM and 
0.125|aM concentrations of ICG, respectively. Fig. 7 plots the lifetime calculation results 
of stage 160 for the two concentrations for DTTCI (hollow symbols) and ICG (solid 

30 symbols) as based on the measurements illustrated in Figs. 5 and 6. It should be 

understood that optical parameters of the samples were not known a priori, and lifetime 
determinations were made without a reference luminophore, being based instead on the 
phase-shift measurements shown. Results of the two-parameter fits of expression (5) to 



70. 
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the phase-shift data are indicated by curves that join the data symbols obtained at the 
excitation and the emission wavelengths. The optical coefficients extracted for the 
samples containing the various concentrations of DTTCI or ICG are summarized in Table 
1. 

5 Experiments were also conducted determining lifetime based on relative 

modulation attenuation measurements with detector 48. Fig. 8 provides relative 
modulation magnitude versus modulation frequency plots of fluorescence, emission, and 
excitation characterizations for the two DTTCI concentrations as illustrated by different 
curve-fitted symbols. Fig. 9 provides relative modulation magnitude versus modulation 

10 frequency plots of fluorescence, emission, and excitation characterizations for the two ICG 
concentrations as illustrated by different curve-fitted symbols. Fig. 10 provides the 
lifetime calculations for the two concentrations of DTTCI and ICG based on the 
measurements plotted in Figs. 8 and 9. 

All publications, patents, and patent applications cited in this specification are 

15 herein incorporated by reference as if each individual publication, patent, or patent 

application were specifically and individually indicated to be incorporated by reference 
and set forth in its entirety herein. Furthermore, the following are hereby expressly 
incorporated by reference: U.S. Patent Number 5,865,754 to Sevick-Muraca et al.; U.S. 
Patent Number 5,818,583 to Sevick-Muraca et al.; pending U.S. Patent Application Serial 

20 Number 09/297,895 to Sevick-Muraca et al. filed on 30 June 1999; pending U.S. Patent 
Application to Sevick-Muraca et al. entitled "Imaging of Light Scattering Tissues with 
Fluorescent Contrast Agents" filed 6 August 1999 as a national stage application of 
International Application Number PCT/US98/023 54 (U.S. patent application serial 
number not yet available); and U.S. Provisional Patent Application Serial Number 

25 60/103,609 to Sevick-Muraca et al. filed on 9 October 1998. While the invention has been 
illustrated and described in detail in the drawings and foregoing description, the same is to 
be considered as illustrative and not restrictive in character, it being understood that only 
the preferred embodiment has been shown and described and that all changes, equivalents, 
and modifications that come within the spirit of the inventions defined by following claims 

30 are desired to be protected. 
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CLAIMS 

What is claimed is: 

1 . A method, comprising: 

interrogating a light scattering medium including an amount of a selected 
luminophore with light at an emission wavelength of the luminophore; 

sensing multiply scattered light at the emission light wavelength in response to said 
interrogating to provide a first optical characterization of the medium; 

exposing the medium to light at an excitation wavelength of the luminophore; 

detecting a multiply scattered luminescence at the emission wavelength in response 
to said exposing; and 

determining a value corresponding to lifetime of the luminophore from the first 
optical characterization and the luminescence. 

2. The method of claim 1 , further comprising illuminating the medium with 
the excitation light wavelength and sensing multiply scattered light at the excitation 
wavelength in response to said illuminating to establish a second optical characterization 
of the medium. 

3. The method of claim 1 or 2, wherein the first optical characterization 
corresponds to photon density of the medium with respect to the emission light 
wavelength. 

4. The method of any of claims 1-3, wherein the emission light wavelength 
and the excitation light wavelength are modulated at a predetermined frequency and said 
sensing and said detecting include determining at least one of a phase difference and a 
modulation attenuation. 

5. The method of any of claims 1-4, wherein said exposing is performed with 
modulated light of the emission light wavelength, said sensing includes detecting the 
multiply scattered light at each of a number of different frequencies of the modulated light, 
and the first optical characterization includes determining a number of values by 
regression with respect to data obtained from said sensing over the different frequencies. 

6. The method of any of claim 1-5, wherein the first optical characterization 
includes an absorption coefficient and an isotropic scattering coefficient of the medium at 
the emission light wavelength. 
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7. The method of any of claims 1-6, wherein the light scattering medium is a 
living tissue and the luminophore is an exogenous probe introduced into the tissue to 
perform a medical diagnosis. 

8. A method, comprising: 

5 exposing a light scattering medium including an amount of a luminophore to a 

number of different light wavelengths selected relative to the luminophore; 

establishing a plurality of optical characteristics of the medium by sensing 
multiply scattered light at each of the different wavelengths; and 

determining a value corresponding to lifetime of the luminophore from the optical 
10 characteristics and a multiply scattered emission at a first one of the light wavelengths 
caused by exposure of the medium to a second one of the light wavelengths. 

9. The method of claim 8, where the first one ef the light wavelengths 
corresponds to an emission wavelength of the luminophore and the second one of the light 
wavelengths corresponds to an excitation wavelength of the luminophore. 

15 10. The method of claim 8 or 9, wherein the optical characteristics correspond 

to absorption and isotropic scattering coefficients of the medium at each of the different 
wavelengths. 

11. The method of any of claims 8-10, wherein the different wavelengths are 
each modulated over time and said exposing and said establishing are performed for each 

20 of a number of different modulation frequencies. 

12. The method of any of claims 8-11, wherein said establishing includes 
determining at least one of a phase difference and a modulation attenuation. 

13. The method of any of claims 8-12, wherein the light scattering medium is a 
living tissue and the luminophore is an exogenous probe introduced into the tissue to 

25 perform a medical diagnosis. 

14. A method, comprising: 

evaluating a light scattering medium including an amount of a luminophore with a 
light interrogation system, the system including a first device optically coupled to a first 
site of the medium and a second device optically coupled to a second site of the medium 
30 spaced apart from the first site; 

subjecting the medium to a first light wavelength selected relative to the 
luminophore; 
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optically coupling the first device to the second site and the second device to the 
first site; 

exposing the medium to the first light wavelength after said optically coupling; and 
determining a value corresponding to lifetime of the luminophore from a first light 
5 output detected in response said subjecting and a second light output detected in response 
another of said exposing. 

15. The method of claim 14, wherein the first device includes a first light 
sensor, the second device includes a second light sensor, and the first light wavelength 
corresponds to at least one of an emission wavelength and an excitation wavelength of the 

10 luminophore. 

16. The method of claim 14 or 15, wherein the first light output and the second 
light output are each of the first light wavelength. 

17. The method of any of claims 14-16, wherein the first light wavelength 
corresponds to an emission wavelength of the luminophore; and said subjecting, said 

15 optically coupling, and said exposing are performed to provide a light scattering 

characterization of the medium at the emission wavelength, the first output light and the 
second output light being of the first light wavelength. 

18. The method of any of claims 14-16, wherein the first light wavelength 
corresponds to an excitation wavelength of the luminophore; and said subjecting, said 

20 optically coupling, and said exposing are performed to provide a light scattering 

characterization of the medium at the excitation wavelength, the first light output and the 
second light output being of the first light wavelength. 

19. The method of claim 18, further comprising repeating said subjecting, said 
optically coupling, and said exposing with a second light wavelength in place of the first 

25 light wavelength, the second light wavelength corresponding to an emission wavelength of 
the luminophore to provide a light scattering characterization of the medium at the 
emission wavelength, the first light output and the second light output being of the second 
light wavelength. 



30 optically coupling, and said exposing; the first light output and the second light output 
being of a second light wavelength corresponding to an emission wavelength of the 
luminophore. 



20. The method of claim 18, further comprising repeating said subjecting, said 
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21. The method of claim 14, wherein the first light wavelength corresponds to 
an excitation wavelength of the luminophore and the first light output and the second light 
output are of a second light wavelength corresponding to an emission wavelength of the 
luminophore. 

5 22. The method of any of claims 14-21, wherein said subjecting and said 

exposing are performed with a modulated light source spaced apart from the first site and 
the second site. 

23. The method of any of claims 14-22, wherein the first light output and the 
second light output correspond to at least one of a phase difference and a modulation 

10 attenuation. 

24. The method of any of claims 14-23, further comprising selectively filtering 
light according to wavelength. 

25. The method of any of claims 1-24, further comprising introducing the 
amount of luminophore into the medium. 

15 26. The method of any of claims 1-25, wherein the luminophore is a 

fluorophore. 

27. The method of any of claims 1-26, wherein said determining is performed 
in accordance with a relationship corresponding to multiple light scattering behavior of the 
medium. 

20 28. The method of claim 27, wherein the relationship is based on a diffusion 

equation model of multiply scattered light in the frequency domain. 

29. An apparatus including an light source apparatus, light detection 
instrumentation, and means for performing the method of any of claims 1-28. 

30. A system, comprising: 

25 light source instrumentation to selectively illuminate a light scattering medium 

including a luminophore; 

detection instrumentation to detect multiply scattered light output from the medium 
in response to illumination by said light source instrumentation; and 

a processor operatively coupled to said detection instrumentation to determine a 
30 first optical characterization of the medium from a first multiply scattered light output of a 
first illumination light wavelength and a second optical characterization of the medium 
from a second multiply scattered light output of a second illumination light wavelength 
different than said first illumination light wavelength, said processor being operable to 
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calculate a value corresponding to lifetime of the luminophore from the first optical 
characterization, the second optical characterization* and a multiply scattered emission cf 
the luminophore from the medium in response to excitation. 

3 1 . The system of claim 30, further comprising a container to receive the 
medium, that is optically coupled to said light source instrumentation and said detection 
instrumentation. 

32. The system of claim 30 or 3 1, wherein said light source instrumentation 
includes a pair of light sources, a first one of the light sources providing light at the first 
illumination light wavelength and a second one of the light source providing light at the 
second illumination light wavelength. 

33. The system of any of claims 30-32, wherein said light source 
instrumentation is operable to provide modulated light of the first illumination light 
wavelength and the second illumination wavelength. 

34. The system of any of claims 30-33, wherein the first illumination light 
wavelength corresponds to an emission wavelength for the luminophore and the second 
illumination light wavelength corresponds to an excitation wavelength for the 
luminophore. 

35. The system of any of claims 30-34, wherein said detection instrumentation 
includes a first light sensor optically coupled to a first position in the medium and a 
second light sensor optically coupled to second position in the medium spaced apart from 
the first position, said detection instrumentation being reconfigurable to optically couple 
said first light sensor to said second position and said second light sensor to said first 
position. 

36. The system of claim 35, further comprising a first optical coupler to 
transmit light from said first position to said first or second sensor, and a second optical 
coupler to transmit light from said second site to said first or second sensor. 

37. The method of any of claims 30-36, wherein said processor calculates said 
value in accordance with a relationship corresponding to multiple light scattering behavior 
of the medium. 

38. The method of claim 37, wherein said relationship is based on a diffusion 
equation model of multiply scattered light in the frequency domain. 
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39. The method of claim 30-38 wherein said detection instrumentation provides 
one or more output signals to said processor corresponding to at least one of a phase 
difference and a modulation attenuation. 

40. A system, comprising: 

5 means for illuminating a light scattering medium including a luminophore; and 

means for characterizing light scattering behavior of the medium for an excitation 
wavelength of the luminophore and an emission wavelength of the luminophore; 

means for determining lifetime of the luminophore from said characterizing means 
and a multiply scattered light emission from the medium at the emission wavelength in 
10 response to illumination by light at the excitation wavelength. 
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